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Abstract. We study Bose-Einstein condensate in the combined presence of 
time modulated optical lattice and harmonic trap in the mean-field approach. 
Through the self-similar method, we show the existence of sinusoidal lattice 
modes in this inhomogeneous system, commensurate with the lattice potential. A 
significant advantage of this system is wide tunability of the parameters through 
chirp management. The combined effect of the interaction, harmonic trap and 
lattice potential leads to the generation of nonlinear resonances, exactly where the 
matter wave changes its direction. When the harmonic trap is switched off, the 
BEC undergoes a nonlinear compression for the static optical lattice potential. 
For better understanding of chirp management and the nature of the sinusoidal 
excitation, we investigate the energy spectrum of the condensate, which clearly 
reveals the generation of nonlinear resonances in the appropriate regime. We 
have also identified a classical dynamical phase transition occurring in the system, 
where loss of superfluidity takes the superfluid phase to an insulating state. 
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1. Introduction 

During the last two decades, the area of Bose-Einstein condensate (BEC) has seen 
significant advances, both in theoretical and experimental fronts [HE]. Response of 
the BEC changes drastically, when the system is loaded in an optical lattice (OL) 
potential, which is an area of active research miElllElIH]- The observed behavior 
of the system can be divided into two parts: static and dynamic cases. In the former 
one, for the case of deep optical lattices, a metal to insulator quantum phase transition 
|9], driven by the quantum fluctuations can occur, whereas, for shallow lattices, 
the same can occur through a classical dynamical phase transition [B Emni ESI, 
driven by modulational instability. The dynamic case is particularly interesting as 
it leads to a number of interesting phenomena in quantum and classical domains. 
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These include, the observation of dynamical localization [HI [15], Anderson metal- 
insulator tunneling in momentum space [S]. gapped modes in degenerate quantum 
gases HZIIIH], realization of nonlinear coherent modes by means of resonant excitations 
munun] etc. Optical lattices can be made time dependent by temporal variations 
of the depth and the width of the lattice |22|. The lattice can be moved at a velocity 
V = {X/2)6iy', where A/2 is the spatial period of the potential and Sv is the stable 
detuning between the two laser beams, which also allows the lattice to accelerate. 
The amplitude of the lattice potential can be tuned through the following relation: 
Vo « 2hV,]i/Eji, riji being the Rabi frequency and En = 2h?TT'^/mX^, the recoil energy. 

A number of recent studies have investigated the condensate dynamics in the 
presence of time modulated optical lattices. It is worth pointing out in this context 
that, Fallani et. at, have observed dynamical instability in one dimensional moving 
optical lattice [52] . Recently, the region of the dynamical instability of the condensate 
with higher order nonlinearities has been explored in a time dependent lattice potential 
[55] . A detailed analysis has been made in [53], for the creation, evolution and 
manipulation of wave packets in the combined potential of optical lattice and a 
harmonic oscillator (HO) trap. Some time back, Madison et. al. [25], studied the 
quantum transport of atoms in the presence of an accelerating lattice. They observed 
resonances, when a small oscillatory trap was added to the system. The number of 
atoms decreases rapidly when the resonances occur. The spatio-temporal dynamics 
of the condensate has also been investigated in the presence of moving optical lattice, 
using direct perturbation and Melnikov-function methods |26| . The existence of sub- 
diffractive solitons and chaos has been shown to exist for a time dependent shallow 
optical lattice [57]. These authors also showed the possibility of suppression of chaos 
and its control through the intensity of the lattice potential. Recently, nonlinear 
resonance phenomenon has been studied in BEC, where the occurrence of resonances 
has been investigated through periodic modulation of the scattering length [58l 1591 [3(1] . 
Unlike the methods opted in these studies, we take a different approach by introducing 
a lattice potential and show that it is possible to generate the nonlinear resonances in 
BEC, through controlled chirp management. 

In this paper, we present a detailed study of the dynamics of BEC in the combined 
presence of a general time dependent optical lattice and a harmonic oscillator trap, for 
both attractive and repulsive atom-atom interactions. The exact sinusoidal solutions 
are naturally chirped, which allows them to be compressed and accelerated in a 
controlled manner. The presence of the trap drastically changes the dynamics of 
the condensate. We show the manifestation of nonlinear resonances, which arise due 
to the combined effect of harmonic trap, interactions and lattice potential. These 
resonances occur, when the matter wave changes its direction. The density of the 
condensate is found to be maximum at this point. In a harmonic trap, the center of 
mass (COM) undergoes periodic oscillation at trap frequency, as a consequence of the 
Kohn theorem, which makes the matter wave change directions periodically. We have 
shown that by engineering the trap parameter, one can have various COM profiles 
that lead to interesting physics. We then consider the BEC in the presence of optical 
lattice, when the harmonic trap is switched off. In case of a static OL potential, 
the BEC undergoes a rapid nonlinear compression, leading to a resonance behavior. 
This phenomenon of rapid nonlinear compression is analogous to the effective pulse 
compression in nonlinear optical fiber |3Ij . For the sake of completeness, we also 
analyze the system in the presence of an expulsive harmonic trap, where BEC spreads 
out. For gaining a deeper insight of these resonances, we investigate the energy 
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spectrum of BEC both analytically and numerically. In the combined presence of 
harmonic trap and OL potential, the condensate shows a resonant increase in energy 
at certain points in the scaled time variable, which supports our observation through 
the density distribution. In case of static OL, a sudden increase of energy mimics the 
resonance behavior of the system. For comparison, we perform a numerical simulation, 
which clearly shows the generation of nonlinear resonances, similar to the analytically 
obtained results. This system also exhibits a classical dynamical superfluid insulator 
transition (DSIT), where superfluidity breaks down and the condensate transits to 
an insulating phase. This phase transition is of the first order and occurs due to 
modulational instability m- We show the stability of the obtained solutions in the 
repulsive domain, using Vakitov-Kolokolov (VK) criterion [32]. 

The paper is organized in the following way: in Sec. 2, we present the theoretical 
model and derive the governing mean-field equation that describes the system. The 
exact solutions, described in Sec. 3, are found in different regions, describing the 
sinusoidal excitations in the presence of both lattice and harmonic trap. In Sec. 4, the 
nature of these obtained solutions and generation of nonlinear resonances in different 
parameter regimes have been analyzed. The generation of nonlinear resonances is 
discussed in Secs. 4.1 and 4.2. The energy spectrum has been investigated in Sec. 5, 
where we find that it supports the generation of the nonlinear resonances. Analytical, 
as well as numerical results, are illustrated respectively, in Secs. 5.1 and 5.2. The 
occurrence of dynamical phase transition in the system is described in Sec. 6. This 
classical phase transition occurs at the point, where the energy expression becomes 
non-analytic. In Sec. 7, we perform stability analysis using Vakhitov-Kolokolov 
criterion and show that the obtained solutions are stable for repulsive atom-atom 
interaction. We draw our conclusion in Sec. 8 and discuss the possible future 
prospects. 

2. Theoretical model 

We consider here a BEC immersed in a time modulated OL potential and highly 
elongated harmonic trap in the longitudinal direction. In the ultra-cold regime, at 
very low temperature, the mean field Gross-Pitaevskii (GP) equation well captures 
the condensate dynamics: 

- '"W) (1) 

Here, m is the mass of the atoms, Vext{r,t) = ytr{x,y) + Vi{z,t) is the external 
trapping potential; Vtr{x, y) and Vi{z, t), respectively, represent the potential along the 
transverse and longitudinal directions. U{t) = jg the strength of the atom- 

atom interaction and v{t) is the time dependent chemical potential. For obtaining a 
cigar-shaped BEC, one applies a strong oscillator trapping potential of frequency uj±, 
Vtr{x,y) = ^mijj'j_{x‘^ + y'^), along the transverse direction. Assuming tight transverse 
confinement, the trial wave function can be factorized: 'k = tjj{z,t)(l^o{x,y), where 
(j)o{x,y) is represented by a Gaussian ansatz (^[33]: 

I 

(t>o{x,y) = C —"“i , (2) 

y 

with a_L = \J xJLx • From now onwards for convenience, we take ijj{x,y) —^ ^|J 
throughout this paper. In the weak coupling regime, one finds that the quasi one 
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dimensional GP equation takes the following form in terms of the dimensionless 
variable [551 15S11?7| : 

(3) 

Here, Vi{z,t) is the trapping potential along the longitudinal direction. The scaled 
two-body interaction is given by, g{t) = where, kj^ is the wave vector. 

The dimensionless time, spatial coordinate and the wave function are, respectively, 
scaled as, t' — >• {hk\/m)t, z' —>■ k^z and ip^z^t) —>■ 'ip{z,t)/V^ [33 (for notational 
convenience, we replace t' by t, z' by z and tp' by ip in Eq. ra). The chemical potential 

v{t) has been scaled in terms of the recoil energy Eji = such that = Efjy(t). 

Keeping in mind the fact that most of the atomic systems are inhomogeneous 
due to the presence of the magnetic or optical traps, we consider the BEG in 
presence of the harmonic oscillator potential, in addition to the optical lattice: 
Vi{z,t) = Vb(t) cos^(T(z, t)) -|- \M{t)z'^, with T(z,t) = k{t){z — l{t)). For notational 
convenience, we use T{z,t) as T, throughout this paper. Here, Ipt) = 
is the position, with v{t) being the COM motion of the condensate and k{t) is the 
inverse of the width. Since, the potential Vi{z,t) is scaled by the recoil energy Er, 
the time dependent amplitude of the lattice potential Vo(t) and frequency of the 
harmonic trap M(t), are both scaled by the recoil energy Eji {Vo{t) = VQ{t)Ei^ 
and M{t) = M'{t)Eji, and henceforth, we replace Ho(t) —>■ Vo{t), M'{t) —)• M{t) 
in the expression of Vi{z,t) above). Control of BEG and its solitonic and sinusoidal 
excitations through the time dependent trap parameter M(t) is an area of signihcant 
current interest |38l [39l |40l El E]- It has been observed that, the center of mass 
of solitons, as well as its width and amplitude get coupled to the trap parameters, 
through which they can be accelerated or compressed [131 E], which is useful for 
coherent atom optics gSlIlSlE]- 

3. Sinusoidal excitations 

A number of methods have been employed in order to solve these type of systems. 
Among them, variational approximation method and perturbation methods are well 
known. In this paper, we follow a different approach known as self-similar method, to 
examine the exact analytical solutions and corresponding dynamics of this system with 
harmonic trap and lattice potentials. The first step of this method is to employ an 
ansatz solution, which consists of a non-trivial phase that is necessarily chirped. This 
type of chirped phase often arises in nonlinear optics and fiber optics as acceleration 
induced inhomogeneity, owing its origin to the interplay between the interaction 
and harmonic trap. We, therefore assume the following ansatz for the BEG profile 

HaBOlE]: 

Pj{z, t) = . (4) 

where, cr(T) is the density and x(T) is the nontrivial phase, controlling the 
supercurrent of the condensate. The chirped phase (p{z,t), is of the form: 
(p{z,t) = —\c{t)z'^. Gurrent conservation yields, where 5 = 

^-\/(WE'2)7+”2a)(T^E^/I). Here, the two-body interaction strength g{t), the lattice 
depth Vo{t) and the chemical potential i/(t) are rescaled by the inverse of the width 
of the condensate k{t), such as, g{t) = nk{t), Vo{t) = ak"^{t) and i/(t) = fj,k‘^{t), 
respectively, k, a and g, are, respectively, the reduced two-body interaction strength, 
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amplitude of the lattice potential and chemical potential. The chirp parameter c{t) 
can be found from the following Riccati equation: 


^-c\t) = M{t). (5) 

The inverse of width k{t) is related to the chirp parameter c(t) by, = k{t)c{t). 

From Kohn’s theorem [48], it is known that for a BEC, confined in a harmonic 
trap with constant frequency, the COM oscillates with the frequency of the trap, 
giving rise to a resonance known as the Kohn mode. It is also known that when the 
trap frequency varies with time, the COM motion is no longer decoupled from the 
trap jlH]. In order to have a complete control over the spatio-temporal dynamics of 
these sinusoidal excitations, we separate the equation governing the Kohn mode gOj: 

+c{t)l{t) = k{t)u. Here, we have introduced a parameter u, which is the velocity 
of the condensate in the absence of the trap. When the harmonic trap is switched off 
{M{t) = 0), the COM moves with a constant velocity u, as is discussed in Sec. 4.4. 
The coupling of the center of mass motion with other degrees of freedom and trap 
parameters, is worth investigating when both trap and lattice co-exist. Interestingly, 
we observe that the above can also be written in terms of the harmonic trap, from 
where the consequences of the Kohn’s theorem can be understood directly: 

^+Mmt)=0. ( 6 ) 

Therefore, it can be clearly seen that the COM solely depends on the harmonic 
trap. Now, the real part of the CP equation, in terms of the density, reduces to 

dasniiiiiisDi: 


QJ^2 O 1 


8\ dT J 
+ acos^(T)a^{T) + 26^ = 0. 


\T) 


(7) 


The first two terms correspond to the contribution from the dispersion term in Eq. 
([^; third, fourth and fifth terms, respectively address to the chemical potential, atom- 
atom interactions and the lattice potential. The last term originates from the current 
conservation and is directly related to the supercurrent of the condensate. For constant 
a, one finds a self-similar solution in terms of the density: 

aiT) = a + bcos^{T). (8) 

where, a = and b = The periodicity of the density modulation is same as 
that of the lattice potential as seen in Eq. Notice that, for a{T) to be positive 
semi-definite, considering the atom-atom interaction to be repulsive, fj, > a + 1/2, 
when a is positive. For, negative a, ji > 1/2. When the interaction is taken to be 
attractive, for positive and negative a, p, > 1/2 and a-\- 1/2, respectively. 


4. Results and discussions 

Resonances, in general, occur when the frequency of the external driving force matches 
with the natural frequency of the system. Resonances are well understood in various 
linear dynamical systems. The investigations of such resonances in nonlinear systems 
are far more exciting, as well as nontrivial. The presence of nonlinearity in the system 
drastically changes the dynamics of such resonances. In this section, we analyze the 
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Figure 1. Fig. (a) represents the density in presence of a regular harmonic trap, 
where at the point of nonlinear compression, the matter wave changes its direction. 
A sharp resonance behavior is observed, when the nonlinear compression takes 
place. Fig.(b) corresponds to the expulsive oscillator, where the BEC spreads. 
Here, the time t and space co-ordinate z have been represented as dimensionless 
parameters after rescaling. In these figures, the parameter values that have been 
used are: cq = 1, fco = O-Si to = 0-4, o = 0.15, u = 0.8, k = 0.8 and g = 1. 


obtained solutions in different parameter regimes and show the generation of nonlinear 
resonances through the time modulation of OL, both in presence and absence of 
harmonic trap, with controlled chirp management. 

4-1. Regular harmonic trap: M{t) = 

In the first case, we assume the trapping parameter as constant: M{t) = which 
represents a regular harmonic oscillator trap. Here, two solutions are possible: 
(i) c{t) = q tan qt and (ii) c(t) = —q cot qt. Consequently, the other parameters 
can be obtained. For (i), k(t) = k^/cos qt and Vo{t) = Vq/cos'^ qt and for (ii), 
k{t) = fco/singf and Vo{t) = Vg/ sin^ qt, with Vg = a/cg. The position of the matter 
wave is oscillatory: (i) l{t) = locosqt and (ii) Ht) = Zgsingf, respectively, with the 
initial position Iq. It is clearly seen that the Kohn mode oscillates with the trap 
frequency q, as predicted by Kohn’s theorem. The condensate eventually evolves into 
periodic oscillations under the modulation of optical lattice and harmonic trap. As 
depicted in Fig. the sinusoidal wave propagates with a periodicity tt and at the 
point of nonlinear compression, the density becomes maximum, wherein, the matter 
wave reverses its direction. Here, the supercurrent takes the form, 

J{t) = k{tmz,t)\^^^=25k\t). (9) 

Fig.§, shows the center of mass profile with respect to time and trap frequency. 
In this case, we consider a general solution of Eq. In this regard, it is worth 
mentioning that one can engineer the trap parameter on demand, in order to make 
the COM profile of a desired form. As an example, we have explicitly checked that 
if the trap parameter varies exponentially {M{t) = e'^*), the COM profile is described 
by Bessel functions of both first and second kind, which are solutions of Eq. ([^. The 
profile is shown in Fig.([^), which shows that the COM no longer oscillates with the 
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trap frequency and shows damped oscillations. One can investigate in detail the effect 
of the trap parameter on COM, which we have left for future investigation. 


0 0.2 0.4 0,6 0.8 1.0 0 0.2 0.4 0.6 0.8 1.0 



Figure 2. COM profile is depicted in Fig.(a), in presence of a regular harmonic 
trap. It shows that the COM oscillates with the trap frequency. Fig.(b) presents 
the COM profile, described by the Bessel functions, when the trap parameter 
varies exponentially. The parameter values used here are same as in Fig. 0 - 


4-2. Static lattice: M{t) = 0 and l{t) =0 

We now consider the case of a static optical lattice potential by putting l{t) =0. For 
effective pulse compression [31] , we consider a time dependent chirp parameter of the 
form, c(t) = and obtain, k{t) = with the initial condition c{t = 0) = cq 

and k{t = 0) = fco- The compression of the matter wave is shown in Fig.(|^). Chirping 
leads to a rapid compression of the BEC proHle and the density tends to a singularity, 
when t —>■ I / cq . It is evident that these self-similar excitations are purely controlled 
through the chirped phase, under certain parametric conditions, which in turn leads 
to efficient compression and amplihcation. The rapid nonlinear compression shown 
here is analogous to the findings by Moores in optical fiber m- 

4-3. Expulsive osdilator:MIt) = —q^ 

For the sake of completeness, we consider M{t) = —which corresponds to an 
expulsive oscillator. The Riccati equation in this case, takes the form: = 

c^{t)—q^. Similar to the previous case, we have two solutions: (i) c(t) = —qtanhqt and 
(a) c(t) = —qcothqt. The inverse of width and lattice amplitude for (i) are: k(t) = 
kg/coshqt, Vo(t) = Vg/ cosh‘s qt and for (ii) k{t) = fco/sinhgt, Vg{t) = Vo/sinh^gt. 
The density profile for this case is shown in Fig. ([^). As time increases, the amplitude 
of the propagating wave decreases, making the BEC spread out. The distortion in the 
density profile, as seen in Fig.([^), is due to the presence of the trap parameter. 

4.4- Moving lattice: M[t) =0 

Here, we discuss the dynamics of BEC in an optical lattice potential, when the 
harmonic trap is switched off, where the condensate has no acceleration and hence 
moves with a constant velocity. The inverse of width and the amplitude of the lattice 
potential have been allowed to vary with time. In this case, the chirp parameter c{t) 
can be determined by taking M(t) = 0 in Eq. = c^(t). In this scenario, 

the solution retains its original form. The position of the matter wave can be found 
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Figure 3. Fig.(a) depicts the sinusoidal BEC profile in presence of an expulsive 
oscillator, where the chirped BEC undergoes expansion. The distortion in the 
density profile appears due to the presence of the trap. Fig.(b) corresponds to a 
parameter domain when harmonic trap is switched off (for moving lattice), where 
BEC undergoes expansion, without any distortion. The parameter values used 
here are same as in Fig. 0 - 


from the relation: l{t) = Iq(1 + cpt), showing a constant velocity v{t) = IqCq = u. 
It is interesting to note that in the absence of harmonic trap, the center of mass 
motion no longer oscillates, but follows a linear profile with time. In the case, with 
c{t = 0) = —XjcQ, we find c{t) = —l/{t + cq ) and k{t) = koCo/{t + cq ). As time 
increases, the inverse of width, as well as the amplitude of the matter wave decreases, 
leading to the spreading of the BEC, as depicted in Fig.([^). In this case, since the 
harmonic trap is switched off, the density profile is smooth and no distortion has been 
seen. 


5. Energy of the excitation 


5.1. Analytical results 


To explore the dynamics of this parametrically driven system, we study the temporal 
behavior of energy, both in the absence and presence of the trap. The energy of the 
condensate can be obtained from the following equation: 


E = 


[ dz 

1 

dtp 

J 

2 

dz 


+ i.Vl{z,t) - 




( 10 ) 


with the terms in the right-hand side, respectively, representing its kinetic energy, the 
potential energy due to both harmonic and optical lattice along with a contribution 
from the chemical potential and interaction energy. The expression for the energy is 
too cumbersome to list and hence, we concentrate on the energy spectrum directly. A 
closer look at the energy spectrum in Fig.Q, reveals that it actually mimics the density 
distributions, shown in Figs.Q. As seen in Fig.(|^), the condensate shows a rapid 
nonlinear resonant increase in energy in the presence of a harmonic trap and lattice 
potential. At certain values of the scaled time variable, it undergoes rapid nonlinear 
compression, which in turn mimics the occurrence of resonances in this system. These 
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Figure 4. Fig. (a) shows the nonlinear resonances in BEC in presence of both 
harmonic oscillator and lattice potential. Inset depicts the full energy profile, 
which shows the occurence of resonances periodically. In Fig. (b), it is shown 
that resonance also occurs when the harmonic trap is switched off. This is akin 
to the nonlinear pulse compression, seen in m- Fig. (c) shows the nonlinear 
resonances obtained from numerical investigation, when both harmonic trap and 
lattice potential are present. Fig. (d) shows the numerical result of the nonlinear 
resonance, when trap is switched off. The energy has been scaled by the recoil 
energy En. The parameter values considered here are same as in Fig. (f). 


nonlinear resonances occur periodically at the point of nonlinear compression of BEC, 
where the density takes its maximum value. The contribution of the quadratic chirped 
phase to the kinetic energy is solely responsible for this phenomenon. In other words, 
the resonances occur when the driving frequency of the potential matches with the 
natural frequency of the system. The matter wave changes its direction at the point 
of nonlinear compression. This is akin to recently observed resonant behavior in an 
optical lattice, where energy transfer takes place between two bands m- It is worth 
mentioning that the nonlinear resonance also takes place when the harmonic trap is 
switched off. Fig.(|^) shows an effective nonlinear compression when the harmonic 
trap is switched off. Rapid increase in energy is found with time for this driven system 
around t ^ 0.8. The thousand-fold increase in energy seen in both Figs.(|^) and (|^), 
is due to large nonlinear compression. This is analogous to the findings in |3Ij . in the 
context of effective pulse compression in nonlinear fiber optics. 


5.2. Numerical results 


Since, the main result of this paper is to show the generation of nonlinear resonances 
in a BEC subject to an external OL potential, in this section we will only concentrate 
on the parameter regimes, where we observe the occurrence of the nonlinear resonance 
analytically. In our numerical simulation, we have used the split-step Crank-Nicholson 
(CN) method in order to solve the time dependent GP equation. In particular, we 
study the time evolution of GP equation and using Eq.(IO), we compute the energy. 
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In this simulation, we use the code realtimeID.F, published in [52j and modify it 
as per the need of the present system. In our case, we choose dz = 0.01 and 
dt = 0.00001, which is well satisfied under the CN scheme: dz'^/dt < 1. We have 
introduced the nonlinearity g{t) in NSTP = 10000 iterations. Due to the sinusoidal 
nature of the solutions obtained analytically, we consider a periodic function as an 
initial wavefunction, as well as periodic boundary conditions. The nonlinearity is then 
slowly introduced into the system in each time iteration until the desired value is 
attained. As a result, this leads to the hnal wavefunction, after NRUN = 900000 
iterations [35] ■ We then use Eq.(10l to compute the energy. To perform a systematic 
numerical study of Eq. ([^ and consequently Eq.(10l, we use the same parameter 
values as provided in the analytical case. The numerical results are shown in Figs.([^) 
and &)■ Fig.(|^) depicts the energy as a function of time, when both harmonic trap 
and lattice potentials are present. It shows a rapid periodic increase in energy at 
certain time interval. This is exactly similar to that we observed in our analytical 
results, shown in Fig.(|^). In the analytical case, as shown in Fig.([^), the resonances 
occur at t ^ 5 and t ~ 7.6, whereas, numerics shows the occurrence of the same 
around t ~ 4.8 and t ~ 7.8. In simulating the case for static OL, we use same NSTP 
to introduce the nonlinearity and consider NRUN = 100000 for the convergence of the 
wavefunction. Fig.(|^) shows the resonance increase in energy, when the harmonic 
trap is switched off. A sudden increase in energy is seen around t ^ 0.9. This is 
also similar to our analytical findings as shown in Fig.([^). In the analytical case, the 
resonance takes place around t ~ 0.8, whereas, numerics shows the same rapid increase 
in energy around t ^ 0.9. In both these cases, we observe a small difference at the point 
of the occurrence of nonlinear resonances, when comparing with our analytical results 
with numerics. This small difference occurs due to the assumption of the periodic 
boundary conditions in our simulations. The analytically obtained solutions through 
the self-similar method, do not follow the periodic boundary conditions. Therefore, we 
conclude that the spatial periodic boundary conditions are responsible for this small 
deviation. However, this does not alter the generation of nonlinear resonances, which 
is the main result of this paper. Therefore, we conclude that the results obtained from 
numerical simulation exactly support our analytical results. 


6. Dynamical superfluid-insulator transition 

We now examine the possibility of dynamical phase transition in this system. In 
case of shallow lattice, a dynamical superfluid-insulator transition (DSIT), driven 
by modulational instability, was predicted by Smerzi et. ah, using the mean-held 
discrete nonlinear Schrodinger equation |T0]. In the following year, Cataliotti et. ah, 
experimentally observed this classical phase transition in the presence of a stationary 
optical lattice m- It has been found in [T3] that, DSIT can also occur when 
both two and three-body interactions are present. Recently, Fallani et. ah, have 
observed this dynamical instability in one dimensional moving optical lattice |22) . 
Since the dynamical phase transition occurs at the point, where the energy of the 
system becomes non-analytic (NA), we have investigated the same from the exact 
energy expression for all the three cases, separately. As mentioned earlier, the energy 
expression is too long to be reported here, we will only concentrate on the non-analytic 
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term, which in all the three cases is same and given by, 

4P(t)(P(t)a(a + b) - 5^) 


y/a(a + b) 


\/a + h 


The energy is non-analytic at the following two points: a + b = 0 and a = 0, 
which respectively, lead to 2/x + 2a — 1 = 0 and 2/r — 1 = 0. Both of these give 
(5 = 0, which reflects from Eq.([^ that the supercurrent is zero. For the density to be 
a finite positive quantity, we exclude the second point. Taking into account the first 
condition, the supercurrent vanishes at this point and the superfluid phase transits to 
an insulating phase, whose wave function is given by. 


= \/^^^sin(A:(f)2)e*^("’*). (12) 

V K 

Since, the superfluid phase is associated with a finite phase factor x(T), the number of 
atoms on each lattice site are unknown. Therefore, the atoms are free to move within 
the lattice and can easily tunnel from one lattice site to another. The phase coherence 
in the superfluid phase can be experimentally observed through the formation of the 
interference pattern. The modulational instability, occurring in the system, leads to 
a loss of phase coherence and hence the superfluid phase transits into the insulating 
phase. The insulating phase is not characterized by a definite phase and hence, the 
number of atoms per lattice site is fixed. Since the phase coherence is lost in the 
insulating phase, no interference pattern will be formed. 


7. Stability analysis: VK criterion 


For a positive semi-definite density of the sinusoidal excitation in NLS-type equations, 
the stability can be examined through the Vakhitov-Kolokolov criterion [52] ■ It states 
that a necessary stability condition is a positive slope in the dependence of the chemical 
potential on the number of atoms per lattice site N] i.e., dN/dfi >0. If ^ < 0, 
the solutions are unstable and are marginally stable, if ^ = 0 [55] |S3]. It is worth 
mentioning that, in some cases, e.g., for gap solitons obtained in [55|, the VK criterion 
does not apply. The stability of the gap soliton families follow an “anti-VK” criterion: 
dfi/dN > 0. Therefore, in order to obtain the stability condition, we calculate the 
number of atoms per lattice site, 


N = 


— a — 1) a 


sin(2/c(t)7r), 


2k Attk 

which can be controlled through chirp management, 
solutions, found in the previous section. 


(13) 


We observe that for all the 


dN _ k{t) 
9/r K 


(14) 


showing that the stability of the obtained sinusoidal solutions depends on the chirp 
controlled inverse of the width of the condensate, as well as on the nature of the atom- 
atom interaction. Therefore, we conclude that the solutions are stable for repulsive 
atom-atom interaction (k > 0). 
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8. Conclusion 

In summary, we present a detailed study of the parametrically forced BEC in the 
presence of a general time modulated optical lattice and a harmonic trap, with both 
attractive and repulsive two-body interaction. The interaction between atoms and the 
lattice leads to sinusoidal excitations, which can be controlled by chirp management. 
We show through the density distribution and energy spectrum that BEC can be made 
to accelerate and undergo rapid nonlinear compression in the presence of a harmonic 
trap. A resonant increase in energy is observed at certain points of time, when the 
system couples to the harmonic trap. The combined effect of lattice and trap leads 
to the generation of these nonlinear resonances at the point, where the matter wave 
changes its direction. For the static case, when the harmonic trap is switched off, BEC 
undergoes a rapid nonlinear compression, analogous to an effective pulse compression 
in nonlinear fiber optics. Both these observations of resonances are well supported 
by numerical simulation. We also show that a classical dynamical phase transition 
occurs in the system, where the loss of superfluidity takes the superfluid phase to an 
insulating phase. The stability of the obtained solutions is investigated using well 
known VK criterion. We find that the solutions are stable for repulsive nonlinearity 
in the present case, though it strongly depends on the nature of the chirped phase. 
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